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Abstract 

The order contribution to the pressure of massless quantum electrodynamics at 
nonzero temperature is determined exphcitly. An identity is also obtained relating a 
gauge-invariant piece of the pressure at order e^"""*"^ > 1) (from diagrams with only 
one fermion loop) to the pressure at order e^". Prospects for higher order calculations 
are discussed and potential applications are mentioned. 
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The three-loop contribution to the equation-of-state of massless quantum electro- 
dynamics (QED) at temperature T was obtained recently thus extending the well 
known two-loop result 0. While three-loop calculations in cold plasmas (T = 0, but 
nonzero chemical potential /j,) were done some time ago 0, same order calculations 
for hot plasmas have been delayed by technical complications due to the presence 
of statistical distribution factors in loop integrals. The three-loop calculation in hot 
QED [0] was stimulated by in part by a ground breaking three-loop calculation in hot 
0^ theory |]. 

The purpose of this letter is to further the results of by determining explicitly 
the next order (e^) contribution to the pressure of QED with N massless electrons 
at nonzero T (but fi = 0). Higher orders will also be discussed qualitatively. The 
conventions here are similar to [Q and are reiterated for clarity : The imaginary 
time formalism is used, whereby the energies are discrete and imaginary, go = inrcT, 
n being an odd (even) integer for fermions (bosons). Ultraviolet singularities are 
regularised by dimensional continuation (4 D) and renormalisation is through 
minimal subtraction. The wave- vector, Q^^ = {qo,q) lives in a space with Minkowski 
metric, g^i, =diag(l, —1, —1.... — 1) and the measure for loop integrals is denoted by 



J {dq} ^TJ2 J (dq) 



go, odd ' 



J[dq]^T E J{dq) 



qo,even 



where 



Spinor traces are normalised according to Tr (7^7^) = 45f^,^, and the photon is in the 
Feynman gauge, D^^{Q) = g^„/Q^. 

Recall that in a QED plasma static electric fields are screened with an inverse 
screening length given to lowest order by eT/ ^/S. Consequently a naive perturbative 
expansion using bare (unscreened) propagators produces severe infrared (IR) diver- 
gences in diagrams such as in Figs. 1,2. When these diagrams are resummed one 
obtains a finite result of order f^, |^ which is nonanalytic in the coupling e^. Just 
as the term is the "plasmon" correction to the two-loop contribution obtained 
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by dressing the photons, the order term to be presently calculated is the plasmon 
correction to the three- loop result. The object that is required in the analysis 
below is the renormalised, static, one-loop photon polarisation tensor Il^^{qo = 0,q). 
Gauge invariance, Q^U^'^^Q) = 0, implies njo(0,g) = (this is of course true to all 
orders) while explicit calculations yield llij{0,q 0) = O(e^g^) and 

noo(0,g) = m2 + eV%), (2) 

where e is the renormalised coupling. The following limits will be required : 
m^{D ^ 4) = and 

/i^-^ 2 iV (D - 2)/2 N 
h(^ = ^^ = s + 6vr^(D-4) - 

The second term in (^) is the UV counterterm and fi is the mass scale of dimensional 
regularisation. The integral /2 in (|^) is defined through 

fn - / 7^ = (2^"^^-'^ - 1) &n , 



bn 



(Q2)n 

[dQ] 

C(2n+ 1 ~D)V 



2 (-1)" T^^2n ^ V ^ f2n + l-D 



(27r)2" V{n) 

Consider now the three-loop, order e^N , diagram Gi shown in Fig. 3a. It is given 

by 

Gi -e^N r {dK}[dQ dP] Tr{^^ ^7,D"-(P)(^- ^^,D'^'^{Q){^- Q)) 



fi^-^D 2 J K\K-QY{K~Py 

(4) 

From the behaviour of the polarisation tensor discussed above, it is deduced that 
when one or more self-energy subdiagrams are inserted along the photon lines of Gi, 
IR divergences occur only for the noo(0, 0) = m? insertions . These may be resummed 
into an effective propagator 

— il^ _|_ 9/^0 9v0 ^qofi 

^ D,,{Q) + D%{Q). (5) 
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When the effective propagator (|) is used for the Q- photon in Gi the correction 
D*{Q) makes the g-integral of order e, / ^2(^2 ~ tti ~ e, and therefore its contri- 
bution to Gi is O(e^). Similarly if both photons are dressed, the extra correction to 
Gi will be of order e^. Thus the order contribution from dressing Gi is obtained 
by keeping one photon bare, using the correction D* for the other, and multiplying 
the result by two : 

4 / ^ rn^ S,,,o f {dK}[dP] Tr{^, ^7a(^- V)^ ^7o(^- Q)) 



-e^NT J {dq) 



li^-'^D J ' ^'q^{q^ + m?) J K^P^{K - Qy{K - P)^ 

Now scale q = mx in the above equation. This results in an overall external factor 
e^m ~ and since the K, P integrals are infrared safe in the m —>■ limit, one 
obtains the exact contribution from 6G1 as 

_ _e4^y^D-3 f (dx) f {dK}[dP] rr(7o ^7.(^- ^)7" ^7o ^) 



(6) 



^8-20 y a;2(a;2 + l) J K^P^{K - Pf 

(7) 

Notice that for the calculation, the original complicated three-loop integral has 
factorised completely into the product of a finite integral and a relatively simple two- 
loop integral. A similar analysis for diagram G2 (Fig.3b) gives for the contribution, 

G25 -e^NTm^-^ r (dx) f {dK}[dP] Tr{^o ^la{^- Vhoi^- ¥) 



/i8-2« 2 J x^{x^ + l)J K^P^{K-PY 

For the sum of and (H) I obtain 

^15 + ^25 = -e^m^-3T/-2^iV5 2(2-D)(6i-/i)(4-D)/2 (9) 

_g5y4^3/2 



0(0-4). (10) 



Pleasantly, while the individual expressions 6*15 and G25 require the evaluation of 
some two-loop integrals, their sum depends only on the trivial one-loop integrals 
bn and /„, defined previously, and S = J ^^l^qrxy- The final result as ^ 4 in 
(p!0|) is finite as required since the contribution of diagrams coming from electron- 



wavefunction and vertex renormalsations cancels as in the calculation because 
of the Ward identity Zi = Z2, and the restriction to massless electrons. 
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The order e'^N'^ diagram is shown in Fig.l. Inserting self-energy subdiagrams 
along the photon lines produces the equivalent set of diagrams shown in Fig. 2. 
Power counting now indicates that only if the self-energies in Fig.2 are noo(0, q) can 
one obtain contributions of lower order than e^. These insertions sum to 



T 
2" 



(dq) 



In (l + H2£iM) ^ _ noo(0,g) ^ 1 



'noo(0,g)~ 



q^ 



q^ 



(11) 



However (|TID contains as a leading contribution part of the lower order piece due to 
noo(OO) = m^. Subtracting from ([Tl|) the same expression but with noo(0, q) replaced 
by gives, upon using (El), 



T 



{dq) 



In 1 + 



e q h{q) 
q2 _|_ ^2 



2q 



(12) 



This contains all terms of order e^N^/'^ but also many subleading terms. Again the 
scaling q = mx helps to identify the exact pieces. Dropping subleading terms one 
finds 

I I rlTA 'r^ hliW \ 

r2 



G 



35 



(dx) x' h{0) 



x2 + 1 



1 ^ 1 



x^ 



x^ 



(13) 



_e5T4^5/2 ^_i + in(4/ 



^7rV27 



127r2 



+ 



g52-4^5/2 in(T//i) 



(14) 



The expression (|T3D has a simple interpretation and could have been written down 
by inspecting the diagrams in Fig. 2 : It is the subleading contribution from the sum 
of infinite diagrams when, for each diagram, the subleading e^g^ piece is kept of 
exactly one self-energy while the leading e^T^ pieces are taken from the remaining 
self-energies. The e^iV^/^ ln(T/ fx) term in the answer ([I^) is a remnant of wavefunction 
renormalisation. This logarithm cancels against a similar term that is generated from 
the lower order e^iV^/^ plasmon term when the pressure is written in terms of the 
temperature dependent coupling 



(15) 



Alternatively one may also eliminate such logarithms by simply choosing fi = T. 

The fine-structure constant at temperature T is a(T) = e2(T)/4'7r. Defining 
g' = a{T)N/7T, the pressure of QED with massless Dirac fermions at nonzero 
temperature, T, then follows from Refs. 0, |I| and eqns.(|TO[), (|1^ and ([T5|), 
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with 



^ = aQ + g^a2 + g^a^ + g\a^ + h/N)+g\a^ + h^/N) + 0{g''), (16) 
«o = ^(1 + 7^^), (17) 



45 ' 4 

57?^ 



a2 



72 ' 

= 97S' '''' 

04 = -0.8216 ±10^\ (20) 

64 = 1.2456 ± 1.5 X 10"\ (21) 

7r2[l-7-ln(4/7r)l , , 

05 = — ^ ^ I = 0.11473... , (22) 



9^3 



- - 2V/3 ■ ^''^ 
Some observations are made on the short perturbation series above. For = 1 
the fourth order (51^) coefficient (04 + 64) = 0.4240 ±0.00025 is positive^ while the fifth 
order coefficient (05 + 65) = —2.734... is negative. Thus there is a sign ahernation 
separately within the even orders {g^, g'^, g^) and the odd orders [g^, g^). On the other 
hand if one considers a double expansion in g and 1/A^ then the (1/A^)° series {o„,} 
shows sign alternation and so does the (1/A^) series at least for the available 

terms. I do not know if this empirical regularity is suggestive of a general result. 

In T = QED the large order behaviour of perturbation theory with respect to 
the expansion parameter was studied many years ago 0. By contrast, at T 7^ 0, 
resumming the IR divergences creates an expansion in and to my knowledge an 
asymptotic analysis in this case is lacking. Some information about the perturbation 
series can be obtained in the large N expansion. To leading order in 1/A^ (the {o„} 
series), the diagrams which contribute are precisely those in Figs. 1,2 (plus the order 
diagram not shown) with full one-loop self-energy insertions. Consider now only 
the odd terms {o2n+i, n = 1,2, ...}, these are the plasmon (e^"'^^) contributions in the 
large N limit. From power counting as before one deduces that all the 02n+i terms 
are contained in the expression 



T f d?q 



2 J (27r)^ 



A ^ noo(0,g)\ noo(0,g) 



q2 I q2 



(24) 



'This value may be very closely approximated by 3^2/10 — 0.42426...; also note that 04 
-7rVl2 = -0.8224... 
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Here noo(0, q) is the finite one-loop self-energy of eq.(0) at D = 4 and renormalisation 
scale /i = T. A systematic expansion of ( p4D in e determines the series {a2n+i} (parts 
of a2n are also contained in the expansion) and it might be possible to determine 
a2n+i in closed form for any n. 

The full calculation was simple because of the factorisation of the diagrams 
Gi and G2. In general, let the gauge-invariant contribution to the pressure at order 
g2n+2 > from diagrams with one-fermion loop, be -P2n+2- Then the order P2n+3 
contribution is obtained by dressing and using the factorisation property : 



piF ^ f idq)m^ 5P^^^^ 

J q\q^ + m^) 6Doo{0,q^O) ' ^ ^ 

„7-/Yl/2 xplF 

= 2±J1 ^^2«+2 y (26) 

n;,r'"^^(0,0), (27) 



„7^2 Tvrl/2 

^-L plF;2n+2 
g32-2jyl/2 g2plF 



2 Im=0 5 



n > 1 . (28) 



IPI 



Equation ( p7|) follows from the earlier equation upon using the relation [H (^^^ 
T 11^'^ /2V , where V =volume, IPI refers to one-particle- irreducible and 11 is the 
full self-energy. The final equation (P5| ) then follows from the identity |^ noo(OO) = 
e^dP^ /dfi^ . The result (^8]) is that mentioned in the abstract. It is valid^ for the case 



of massless electrons at zero chemical potential fi. For the case n = 1 one may 
verify that ( |28D is satisfied by computing the right-hand-side from known two-loop 
results p, and comparing it with the left-hand-side given by (p!0|). 

Eg . (|28|) suggests that the contribution should be calculable because knowledge 
of three-loop diagrams is at hand |I[ (diagrams with more than one fermion loop 
are not covered by ( pHD but they are easier and are handled as for the case Fig.l 
discussed above). If extra cancellations occur, as in the computation, then the 

calculation might reduce to a two-loop calculation. Note that the factorisation 
is successful because fermions in imaginary time have IR safe propagators. If one 
attempted to calculate the contribution to the pressure of A^0^ theory, the dia- 
gram in Fig. 4 is involved. Now the exact factorisation fails because the two-loop 

•tAn alternative derivation using the ring-summation formula (Q), and the extension to massive 
electrons at nonzero chemical potential will be presented elsewhere. 
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subdiagram has a logarithmic singularity coming from the zero modes of the bosonic 
propagators. Nevertheless the calculation in this case is still feasible because only 
the leading T^(lnfc+ constant) piece of the two-loop self-energy subdiagram (with 
inflowing momentum {0,k)) is needed to determine the order T^A^(lnA-|- constant) 
contribution . 

Calculations similar to those discussed in and here have been performed to 
obtain directly the screening masses in hot QED to high-order 0, and future appli- 
cations to quantum chromodynamics are envisaged. One speculative application 
for QED will be mentioned : there has been much discussion in the literature (see, e.g., 
||Tl| and references therein) concerning a possible strong-coupling phase (at T = 0) 
of QED with an ultraviolet-stable fixed point. Usual perturbative arguments indi- 
cate that at exponentially high temperatures the QED coupling will be strong. If a 
temperature driven transition to a stable nonperturbative phase is possible, then a 
resummed high-order perturbation series [0] for hot QED might be useful for some 
studies. 

In conclusion, some of the questions raised in the last three paragraphs are left 



for future investigations |12 
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Figure Captions 

Fig.l : 

The order e^N'^ diagram. It has an infrared singularity which contributes to the lower 
(e^) order. In Figs. 1-3 the wavy line is the photon. 
Fig.2 : 

Diagrams generated from Fig.l by photon-polarisation insertions. 
Fig.3 : 

The order e^N diagrams which produce e^N^^"^ terms when the photons are dressed. 
Fig.4 : 

The nontrivial three- loop diagram of A^0^ theory studied in |Q. It will contribute a 
In A piece to the free energy when dressed. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9406318vl 



